Seeding example for 1D environment of S&H 2012 #1396
Conversation
| extensive_attributes, | ||
| seeded_particle_index, | ||
| seeded_particle_multiplicity, | ||
| seeded_particle_cell_id, | ||
| seeded_particle_cell_origin, |
There was a problem hiding this comment.
what is the cell origin and position in cell? where are these used?
There was a problem hiding this comment.
These are used to define and track positions of particles in the domain (only used for >0D environments, so do not apply in the case of parcel env)
There was a problem hiding this comment.
I will use the particulator.environment argument to call the seeding function with default arguments when seeding in the parcel environment
| ) # TODO #1387: does not have to be the same? | ||
| v_dry = settings.formulae.trivia.volume(radius=r_dry) | ||
| self.seeded_particle_extensive_attributes = { | ||
| "water mass": np.array([0.0001 * si.ng] * settings.n_sd_seeding), |
There was a problem hiding this comment.
should we use equilibrate_wet_radii rather than prescribing a water_mass for the seeds?
That is due to the However, this made me wonder: why are the number of SDs decreasing with time initially even in the no seeding case? |
This is a good question... I thought it was just precip, but the sd_count is decreasing even before precip onset (and actually faster in the initial time than during the precip period). Is there some spin-up that's happening? I don't know why that wouldn't conserve SDs though. |
I expect it to be caused by the updraft which is only present during the first |
Oh, this totally makes sense. Thanks, Sylwester! |
Update toward solving a subset of goals raised in #1387.
Specifically, I included functionality for selecting
cell id,cell origin, andposition in cellof the seeded super-droplets (SDs) in Shipway and Hill's (SH2012) kinematic driver environment. For this I had to modifyseeding.py,particulator.py, as well as the backendseeding_methods.pyfiles in a way that does not seem to be straightforwardly backward compatible with the parcel model example that @slayoo introduced earlier.I also made another major assumption: rather than simply using a default value of
water massfor the seeded SDs, I update it using thevolumeattribute derived from the wet, rather than dry, radius of seed SDs at the time of injection. If I use the dry radius, the resulting precipitation is significantly lower. This might be a subtle point that is worth discussing in some detail, even if it is for my own personal clarification.Here's the .ipynb file containing the full example.